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We test the scaling behaviour of Wilson, hypercube, maximally twisted mass and overlap fermion 
actions in dynamical simulations of the 2-dimensional massive Schwinger model. We also present 
possibilities to simulate overlap fermions dynamically by replacing the exact overlap operator by 
an approximate version. This approximation is used either as only the guidance Hamiltonian, 
keeping the exact overlap operator as the accept/reject Hamiltonian or for both, the guidance and 
accept/reject Hamiltonian in the Hybrid Monte Carlo algorithm. In the latter case we test whether 
the approximation error can be corrected for by computing the determinant ratio of the exact and 
the approximate overlap operators stochastically. 
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1. Introduction 

Determining the scaling properties of physical quantities computed with different lattice ac- 
tions as used presently for dynamical fermion simulations is certainly one of the essential steps to 
perform a controlled continuum limit. Obviously, performing such scaling tests for lattice QCD 
in four dimensions is (at least) very demanding. Therefore, we have chosen the 2-dimensional, 
massive Schwinger model [p as a test laboratory to address the scaling properties of a number of 
fermion actions for Nf = 2 flavours of dynamical fermions. In particular, we will compare standard 
Wilson [0], Wilson twisted mass at maximal twist, hypercube []|] and overlap fermions [^] in 
their approach to the continuum limit. 

In order to study the scaling behaviour, we will fix the scaling variable z = (nif^/p) 2 / 3 to 
z = 0.2,0.4,0.8 where m/ is the fermion mass in lattice units derived from the PCAC relation 
and j8 is the coupling multiplying the Wilson plaquette gauge action used throughout this work. 
Denoting by e the physical gauge coupling and by a the lattice spacing, j3 = l/a 2 e 2 . At each 
of the values of z given above, we compute the pseudo scalar mass and follow its behaviour as 
a function of (decreasing) values of the lattice spacing. Performing finally a continuum limit of 
our results allows us to compare to analytical predictions that are available from approximations 
of the massive Schwinger model which cannot be solved exactly. For a discussion of the scalar 
condensate we refer to ref. [0], where also a more detailed account of our work can be found. 

Another important aspect discussed in this contribution is a description of several ways to 
perform overlap simulations dynamically. We follow the general idea to replace the exact overlap 
operator D ov by some approximate version Dov Prox having an infra-red safe kernel. To correct for 
such an approximation we test to reweight by the determinant ratio det[D ov ]/det[Do Pprox ]. We also 
have tested the possibility to use Do Pprox as the guidance Hamiltonian in the molecular dynamics 
part of the Hybrid Monte Carlo (HMC) algorithm while keeping the exact overlap operator D ov as 
the accept/reject Hamiltonian. 

2. Scaling results 

The pion mass m n and the PCAC fermion mass rapcAC have been computed by standard meth- 
ods, using the exponential decay of the 2-point function of the pseudo scalar operator P and com- 
puting the PCAC relation, respectively. In table [j] we list the techniques we have used for the 
simulations and to compute the correlators numerically for the different fermion actions. We just 
note that for the twisted mass (tm) fermion case we worked at full twist, realized by the vanishing 
of the PCAC fermion mass calculated from pure Wilson fermions. In this case, the operators used 
for the correlation functions need to be appropriately rotated [Q], of course. We refer to ref. [|7|] 
for the simulation details such as the statistics of our runs and the lattice size and for the precise 
definitions of the different actions and the notation we have employed. In this reference we also 
discuss the finite size corrections we have performed for the pion masses when necessary. 

In fig. [T] we show one example of the continuum extrapolation of m n ^J~$ for the four different 
actions we have used. We can see that for all the actions there is a nice scaling with 1//3, i.e. 
with a 2 . Another observation is that for all four actions the continuum extrapolated values are 
consistent, in agreement with universality. Given this finding, we decided to perform constraint 
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Wilson / tm / hypercube 


overlap 


simulations 


standard HMC 


determinant reweighting 


operators 


conjugate gradient solver 


eigenvalues and eigenvectors 
conjugate gradient solver 


currents 


local / conserved 


local 


pseudo scalar correlator 


(PP) 


(PP - SS> 



Table 1: The simulation techniques and the methods to compute the observables for the different actions we 
have used are listed. 




Figure 1: Scaling test of m n y/j3 as function of 1//3 a 2 at fixed value of z = 0.4. 



z=(m /v ^) 2 / 3 
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0.2 


0.4016 


0.43266 


0.406(4) 


0.4 


0.8032 


0.86532 


0.825(5) 


0.8 


1.6064 


1.73064 


1.755(4) 



Table 2: Continuum values for m n yf]i at fixed value of z = (mf^Jfi) 2 ^. We compare our continuum 
extrapolated non-perturbatively obtained values, denoted with NP, with theoretical predictions using a Sine- 
Gordon approach (SG, [K|) and using a large mass approximation (LG, [PR). 



continuum extrapolations using fits linear in 1//3 with a common continuum value for all actions 
at fixed value of z. These final results from our analysis together with theoretical predictions are 
given in table ^[ The theoretical calculations are performed using the Sine-Gordon model [jg], 
denoted as SG in the table, and from a large mass expansion []|] which we denote as LG. Our own 
non-perturbatively obtained values are denoted as NP. As one can see from table ^, with the high 
precision we could reach with our simulations, deviations from these theoretical expectations exist. 
Only for a value of z = 0.2 there is a consistency with the prediction of ref. [Q] while the results of 
ref. [||] may start to describe the data only for z = 0.8 or larger. 
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2.1 Simulations 

For the Wilson, the hypercube and the twisted mass fermion action a standard Hybrid Monte 
Carlo algorithm (HMC) [[K]] has been implemented for our dynamical simulations. For overlap 
fermions, we generated configurations with the pure gauge action and then reweighted all observ- 
ables with the overlap fermion determinant. 

The reason for this procedure is that direct simulations of overlap fermions using the HMC 
algorithm are very difficult and face the generic problem that the kernel operator used in the sign 
function of the overlap operator can develop very small eigenvalues. We have tried to avoid this 
problem by replacing D ov by some operator Dov Prox that is a good approximation to D ov but which 
is safe against these low-lying eigenvalues. 

The general idea is to write 

detD ov = detD^ prox -det 

One way to implement this idea is to use solely Do V prox for the HMC algorithm of the gauge fields. 
The remaining correction determinant ratio /?det could be implemented as a reweighting factor 
in the computation of a given observable. Another possibility is to use Do V pr ° x as the guidance 
Hamiltonian in the molecular dynamics part of the HMC algorithm only while the exact overlap 
operator D ov is used to evaluate the Hamiltonian for the accept/reject step. 

A crucial question in the first of these approaches is, whether an operator Do Pprox can be found 
such that the fluctuations in /?det are small enough to obtain statistically significant results. We 
decided therefore to test this idea by computing /?det on a number of gauge field configurations 
generated in the pure gauge theory atj8=3onaL=16 lattice. 

As approximations to the overlap operator we used the hypercube operator Dh yp , an overlap 
operator with a modified Chebyshev polynomial to evaluate the sign-function and an explicitly 
infra-red regulated operator, where 

+ , (2.2) 

with Dq representing the kernel operator to construct the overlap operator. Out of these set of 
operators, it turned out that only the operator given in eq. (^2f) led to an acceptable size of the 
fluctuations for the determinant ratio 7?det- However, even in this case, when the value of 8 is 
chosen to be comparable to values of the low-lying eigenvalues of the kernel of the overlap operator, 
the fluctuations in 7?det became so large that this operator does not seem to be practical when the 
reweighting technique is used, see ref. [Q] for a more detailed discussion. 

For the hypercube operator it turned out that although the mass parameter was tuned in a wide 
range, 7?det showed always large fluctuations such that no stable result for a physical observable 
could be obtained. There the square root of the kernel is computed by a Chebyshev polynomial P n e 
of the degree n in the interval [£,1]. In the studies of the modified Chebyshev polynomial P^g we 
found that it is possible to decrease the degree n of the polynomial to n < n, while keeping £ = £ , 
by about a factor of three. Although this certainly accelerates the simulations, we also found that 
already a slight increase of the lower bound £ of the polynomial approximation to £ > £ led again 
to large fluctuations of 7?det- 



Do, 

^approx 
'--'ov 



detDr r ° x -#det 



(2.1) 
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Figure 2: The acceptance rate P acc using the hypercube operator as the guidance Hamiltonian in the Hybrid 
Monte Carlo algorithm. P acc is plotted as a function of the bare overlap fermion mass m^ yp for three lattice 
volumes. 



As a second attempt, we used an approximate overlap operator for the molecular dynamics 
part of the HMC algorithm. The setup is then to use D ov for the accept/reject Hamiltonian while 
Dov P1 ° x will be used only for the guidance Hamiltonian. We started to investigate this idea by using 
the hypercube operator as the guidance Hamiltonian. In fig. || we show the acceptance rate as a 
function of the bare hypercube fermion mass parameter mj v on three lattice volumes. The results 
are very surprising. From our previous experience with computing R^t we expected a quite bad 
behaviour of the acceptance rate P acc . In contrast, we find that P acc is rather large, at least for the 
8 2 lattice. Even on our largest lattice (32 2 ) the acceptance rate does not go to zero but assumes 
reasonable values. Fig. || also demonstrates that P acc develops a maximum at a certain value of 
the bare hypercube fermion mass that does not vary strongly with the lattice size. We remark 
that for a certain threshold step size we could not increase the acceptance rate by decreasing the 
step size keeping the trajectory length fixed. This means that at this point we are sampling only 
configurations that are distributed according to the hypercube action. 



3. Conclusions 

In this contribution we have tested four different lattice fermions in their approach to the 
continuum limit in the 2-dimensional massive Schwinger model with Nf = 2 flavour of dynamical 
fermions. At fixed scaling variable z = (mf^J~fi) 2 /' i = 0.2,0.4,0.8 we have computed the pseudo 
scalar mass m K yf]i for various values of j8 = 1 /e 2 a 2 . 

For all kind of fermions used, Wilson, hypercube, maximally twisted mass and overlap fermions, 
the scaling behaviour of m n ^J~$ appears to be linear in a 2 and lead to consistent continuum values, 
demonstrating nicely universality of the results. The precision of our simulation results reveal that 
the analytical calculations do not provide a satisfactory description of the data for all values of z. 



239/5 



Scaling in the Schwinger Model 



Beatrix Pollakowski 



We also presented possibilities to simulate overlap fermions dynamically by replacing the 
exact overlap operator by an approximate version which is used either as only the guidance Hamil- 
tonian or as the full Hamiltonian in the Hybrid Monte Carlo algorithm. In the latter case we tested 
whether the approximation error can be corrected by computing the determinant ratio of the exact 
and approximate overlap operators stochastically. This way turned out to be not successful since 
this (correcting) determinant ratio showed large and un-acceptable fluctuations. 

Changing the guidance Hamiltonian showed promising results, however. Using the hypercube 
operator, which led to large fluctuations in evaluating the determinant ratio R^ t , as the guidance 
Hamiltonian we find reasonable acceptance rates. Unfortunately, P acc drops significantly when 
the lattice volume is increased. If we think of simulations in four dimensions, we therefore have 
to conclude that the hypercube operator will in this case not be a good choice for the guidance 
Hamiltonian. 

Nevertheless, given the fact that the other approximate overlap operators, modified Chebyshev 



polynomial and the operator in eq. (|2.2|), led to much more stable results for the fluctuations of Rdeu 
we believe that using these operators as the guidance Hamiltonian should lead to much improved 
acceptance rates. We are presently testing this idea. 
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